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Abstract Skeletochronology is a method commonly used for estimating the age of amphibians and reptiles in the wild. 
However, the number of lines of arrested growth (LAGs) does not necessarily reflect age in some species. We validated 
the applicability of this method to an endangered eublepharid gecko, Goniurosaurus kuroiwae, then inferred its longevity 
and age structures in wild populations. We classified young geckos into three groups using previously published data 
for early growth: Group | contained hatchlings before the first winter, Group 2 contained hatchlings after the first win- 
ter, and Group 3 included yearlings after the second winter. LAG numbers in these groups were then compared. All 
individuals in Group | possessed a single LAG, which was considered as a hatching line. Most individuals in Groups 2 
and 3 possessed one and two additional LAGs, respectively (LAGI and LAG2), corroborating the notion that LAGs are 
formed annually. A few geckos exhibited fewer LAGs than expected. Analysis of variations in LAG and marrow cavity 
diameter demonstrated that in animals with fewer LAGs, endosteal resorption or fusion of hatching line and LAG] 
had occurred. LAG2 was never lost by endosteal resorption and was identifiable by its diameter. Thus, the age of adult 
geckos could be determined by counting LAGs outward from LAG2. Application of this method to wild populations re- 
vealed that the longevity of this species is not less than 83 months, but that almost all individuals in fragmented habitats 
die before 50 months, suggesting lower population sustainability in such habitats. 


Keywords age structure, endosteal resorption, hatching line, longevity, rapprochement, Ryukyu Archipelago 


1. Introduction program based on various ecological parameters is 
required, along with the preservation of habitats. Although 
The genus Goniurosaurus is a group of ground-dwelling 
eublepharid geckos distributed in the central part of 


the Ryukyu Archipelago of Japan, southern China, 


ecological traits of G. kuroiwae have been relatively well 
studied among reptiles in the same area (e.g., Tanaka and 
Nishihira, 1987, 1989; Dial and Grismer, 1992; Werner 
et al., 2005, 2006), knowledge of the age of adult geckos 
is still lacking. Age determination of adult geckos would 


and northern Vietnam. Goniurosaurus kuroiwae is a 
relict endemic to 11 islands of the Central Ryukyus, 
and is at risk of extinction mostly due to habitat loss 
by anthropogenic disturbance, predation by introduced 
carnivorous mammals, and illegal hunting for pet keepers 
(Ota, 2010). The necessity of conservation of this species 
has therefore been emphasized (Tanaka, 1996; Ota, 
2000; Honda, 2002). For conservation, a species-specific 


enable estimation of the longevity and age structure of 
a population, facilitating the evaluation of population 
sustainability (Ralls et al., 2002; Di Minin and Griffiths, 
2011). 

Skeletochronology, in which the number of lines of 
arrested growth (LAGs) formed in bone is counted to 


7 - , , provide an indication of age, is an efficient and effective 
Corresponding author: Takaki KURITA, a Ph.D. student in Graduate 


School of Engineering and Science, University of the Ryukyus, 
Nishihara Okinawa, Japan, with his focusing on taxonomy, population 
genetics and natural history of reptiles. 

E-mail: mofumofu_monticola@hotmail.co.jp 

Received: 7 July 2013 Accepted: 10 October 2013 


method for age determination of wild animals. It has 
been applied to various taxa of amphibians (e.g., Measey 
and Wilkinson, 1998; Khonsue et al., 2000; Matsuki 
and Matsui, 2009; Liao and Lu, 2011) and reptiles 
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(e.g., Waye and Gregory, 1998; El Mouden et al., 1999; 
Coles et al., 2001; Horner and Padian, 2004). Despite 
its wide application, the accuracy of this method for age 
determination to support empirical research has only been 
validated in a few studies (Matsuki and Matsui, 2009). 
This lack of validation may increase the risk of erroneous 
age estimations (Gibbons, 1976), so the corroboration of 
this method in the species or species group in question is 
important. 

Skeletochronological age determination has another 
matter that we need to take bones from the specimens or 
live animals, and even in the latter case, it causes more or 
less damage to animals. Although the mark and recapture 
method has often been used for this purpose, it is time and 
labor intensive to mark and trace animals, and sometimes 
it is too hard to investigate multiple populations 
simultaneously. Considering this, skeletochronological 
age determination using digit bones is a more realistic 
way for practical use. Although toe clipping is considered 
to have various effects on amphibians and reptiles (see 
Perry et al., 2011 for related review), use of digit bones 
appears to not have any serious long-term effect on the 
fitness on locomotion of lizards and geckos and would 
thus be the best method to use for G. kuroiwae (Dodd, 
1993; Paulissen and Meyer, 2000). Therefore, we 
considered that skeletochronology using digit bones is 
on the whole a cost-benefit balanced way to determine 
age of G. kuroiwae at the minimum expense of animal 
damage, even in terms of conservation. 

The primary purpose of this study is to determine the 
effectiveness of the skeletochronological method for 
age determination of G. kuroiwae. We also attempted 
to establish a reliable way to count LAGs and to clarify 
the relationship between the resulting numbers and the 
age of the animal. The validity of this method was well 
established (see Discussion); therefore, we applied it to 
the wild populations of G. kuroiwae in Okinawayjima to 
clarify longevity, age structure, and individual growth 
based on the estimated ages of the samples. 


2. Materials and Methods 


2.1 Study animal and sampling Tanaka and Nishihira 
(1987), in their surveys on the annual activity of 
G. kuroiwae in the southern part of Okinawajima, 
demonstrated that: G. kuroiwae reaches the peak of its 
seasonal activity in the summer months of June and July, 
and becomes mostly inactive in winter from December 
to February; hatchlings start appearing from September; 
both males and females attain sexual maturity at around 


20 months of age (during spring after the second winter), 
at which their snout—vent lengths reach to or exceed 
75 mm. Gravid females, as recognized by abdominal 
palpation, appear from May to August, and each female 
produces two eggs per clutch. Females have more than 
one clutch in a single reproductive season, with some 
females possibly producing 3—4 clutches (Tanaka and 
Nishihira, 1989). 

Goniurosaurus kuroiwae were caught in the south 
to central parts of Okinawajima, including the Motobu 
Peninsula (Figure 1), during the period from April 2010 
to March 2012. Snout—vent length (SVL) of geckos 
was measured to the nearest 1 mm using a ruler. Sex 
was determined for presumably mature geckos whose 
SVL is > 75 mm (Tanaka and Nishihira, 1987), on the 
basis of the presence or absence of hemipenal bulges. 
The second digit of the left forelimb was clipped from 
each individual, and preserved in 10 % formalin until 
skeletochronological examination. After toe clipping, the 
animals were released to the original sites of capture. All 
field experiments were carried out with the permission 
of the Cultural Property Section of Okinawa Prefectural 
Board of Education be-cause handling of this species is 
strictly regulated by law. 


Figure 1 Map of the insular region of East Asia and enlargement 
of Okinawajima Island (box). Solid and open circles indicate locali- 
ties of individuals collected from the Motobu Peninsula and main 
body of Okinawajima, respectively. 


2.2 Skeletochronology Scales and skin were removed 
from the preserved digit. The digit was soaked in tap 
water overnight, decalcified with 5% nitric acid for 90— 
120 min, and soaked overnight again. The medial phalanx 
was separated from the other two phalanges and sliced 
into cross-sections 20 um in thickness using a cryostat 
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LEICA CM1500 (LEICA, Wetzlar, Germany). The 
cross-sections were stained with Mayer’s hematoxylin 
for 60 min, soaked in water for 30 min, and mounted 
on glycerin. Two fine-resolution cross-sections of the 
phalanx from each individual were chosen for LAG 
observation. The number of LAGs was counted under a 
light microscope (Eclipse E200; Nikon, Tokyo, Japan), 
and images were documented as digital photographs 
using a digital camera (COOLPIX 4500; Nikon). 

A simple and effective method to demonstrate that 
the number of LAGs accurately reflects age is to utilize 
individuals whose exact age had already been determined 
using another method, such as mark and recapture or 
raising in captivity. In this study, no marked or captive 
individuals were used. However, the age of most 
immature individuals could be unambiguously determined 
on the basis of collection date and SVL at the time of 
collection, according to the data on early growth of G. 
kuroiwae provided by Tanaka and Nishihira (1989). We 
therefore divided our samples into the following four age 
groups for comparison of the number of LAGs: Group 
1 contained hatchlings collected before the first winter 
(from September to November) with SVL < 60 mm; 
Group 2 included those collected after the first winter but 
before the first summer (February to August) with SVL < 
70 mm; Group 3 comprised yearlings collected after the 
second winter (February to June) with SVL 70-74 mm; 
and Group 4 contained age-undetermined adults with 
SVL = 75 mm. This gecko hardly appears at all in the 
winter months (Tanaka and Nishihira, 1987); therefore, no 
individuals were collected in December and January. We 
did not assign juveniles with SVL 70-74 mm collected 
in July to November to any of the 4 groups because their 
status as hatchlings or yearlings were ambiguous based 
on collection date and SVL. 

LAGs are generally formed when the growth of the 
animal slows or stops mainly due to external conditions, 
such as temperature, humidity, and/or food availability. 
Because the activity of G. kuroiwae declines in winter 
(Tanaka and Nishihira, 1987), LAGs in this species, if 
any, are most likely to be formed in that season. However, 
in some amphibians and reptiles distributed in tropical to 
subtropical regions, LAGs are formed in the dry season or 
arid summer (Khonsue et al., 2000; Olgun et al., 2005). 
Therefore, we attempted to confirm the LAG formation 
in winter by comparing the number of LAGs between 
individuals collected before and after the first winter 
(Groups 1 and 2, respectively). Then, we elucidated 
the presence or absence of LAGs formed in summer by 
comparing Groups 2 and 3. If no LAGs were formed in 
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summer, the difference in the number of LAGs between 
Groups 2 and 3 would be consistent with the number of 
LAGs formed in a winter, which was determined by the 
comparison between Groups | and 2. 

Even if a single LAG is formed annually, a discrepancy 
between the observed number of LAGs and the actual age 
of the individual may be caused by endosteal resorption, 
which is expansion of the marrow cavity resulting in 
removal of one or more LAGs from the medial side of the 
periosteal layer of the bone. To infer the number of LAGs 
lost by endosteal resorption, we measured the horizontal 
diameters (direction perpendicular to the palmar—dorsal 
direction) of the marrow cavity, resorption line, which 
represents the limit of the extent of resorption, and LAGs. 
We determined direction of cross-sections based on 
their shape and defined horizontal diameter as distance 
between the left and right margins at a half depth of 
second-phalanx cross-section. These measurements were 
made on digital photos of cross-sections, with distance 
between two points on the cross-section, which overlaid 
ocular micrometer, as a reference. The software ImageJ 
1.43u (National Institutes of Health, Bethesda, MD, USA) 
and method of Abramoff et al. (2004) were used in this 
procedure. 


2.3 Application to the Okinawajima populations 
We inferred the age structure of G. kuroiwae based 
on individuals collected from the Motobu Peninsula 
(Figure 1). Continuous forest suitable to G. kuroiwae 
dominates the limestone area in this peninsula, and the 
series of individuals collected for this study can therefore 
be regarded as a sample from a single population. In 
contrast, forests in the south to central areas of the 
main body of Okinawajima are heavily fragmented and 
shrunk. Therefore, individuals from this region cannot 
be considered representatives of a single population. 
However, we pooled these individuals into a single 
sample due to the small sample size. 

We estimated the individual growth of G. kuroiwae 
by applying the von Bertalanffy growth formula, Lt = 
Lmax{1 — exp[—K(¢t — f))]}, where ¢ is age in months, 
Lt is SVL at the time sampled, Lmax is the asymptotic 
average of SVL, K is the growth coefficient, and to 
is the artifact of the model that represents age when 
the average SVL is zero (Akamine, 2009). For this 
calculation, we lumped the samples from the Motobu 
Peninsula and other areas to combine data for all 
collected individuals because no significant difference 
in the von Bertalanffy growth curve was observed 
between the samples from Motobu and the other areas 
in either sex (male, F;,, = 0.04, P = 0.99; female, F; go 
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= 0.14, P = 0.94; see below). The von Bertalanffy 
growth parameters for male and female were estimated 
separately. Juveniles were used for estimations of both 
sexes following Leclair and Leclair (2011). Difference 
in the von Bertalanffy growth curve between sexes was 
tested using the F-test, following Akamine (2009). 


3. Results 


3.1 Characterization of LAGs, resorption line, and 
marrow cavity For the validation study, we used 133 
individuals of G. kuroiwae, including 35 juveniles, 38 
males, and 60 females with SVL ranging from 37 mm 
to 96 mm. In total, 129 individuals were assigned to 
predefined sample groups, Group 1—4 (n = 6, 19, 6, and 
98, respectively). The remaining four were not assigned 
to any of the groups because it was unclear that they were 
hatchling or yearling (see above). 

The shape of the diaphysis of the second phalanx was 
square and slightly concave on the palmar and dorsal 
sides in hatchlings (Figure 2 A). As individuals grew, 
periosteal bone tissue was deposited mainly on the lateral 
sides of the bone and less so on the palmar and dorsal 
sides (Figure 2 B-E). The cross-sectional shape of the 
second phalanx became trapezial with slightly concave, 
longer palmar margins and convex, shorter dorsal margins 
as the animals aged. In larger individuals, the bone grew 
further in the lateral direction and became almost oval 
in shape (Figure 2 F). LAGs were clearly observed as 
dark hematoxylinophilic lines. All but one individual had 
at least one LAG and many elliptical osteocytes in the 
periosteal bone. Especially in adult geckos, a resorption 
line, which was distinguishable from ordinary LAGs 
due to its distorted shape and location between the 
periosteal and endosteal bone elements, was frequently 
observed. In the endosteal bone, which constitutes bone 
tissue medial to the resorption line, some faint lines and 
a few osteocytes were observed, but these lines were not 
counted as LAGs. 

One LAG was evident in all hatchlings collected 
before the first winter (Group 1), which was considered as 
a hatching line. The horizontal diameters of the hatching 
line and absorbed zone (the larger one of resorption line 
or marrow cavity in a cross-section) of animals in Group 
1 ranged from 165.0 um to 207.1 um and from 128.6 um 
to 181.6 um, respectively. 

In hatchlings collected after the first winter (Group 
2), two LAGs were usually evident (n = 15), but animals 
with no LAG ( = 1) or only one LAG (n = 3) were also 
observed (Figure 2 B-D). The diameters of the outer 


and inner LAGs, and absorbed zone of individuals with 
two LAGs ranged from 179.9 um to 253.2 um, from 
159.2 um to 216.3 um, and from 137.3 um to 201.0 um, 
respectively. The diameters of the LAG and absorbed 
zone of individuals with one LAG ranged from 168.1 um 
to 219.4 um and from 149.5 um to 183.4 um, respectively. 
The diameter of the absorbed zone of the individual 
with no LAG was 145.5 um and 154.5 um in two cross- 
sections. 

Juveniles collected after the second winter (Group 3) 
had two (n = 1) or three LAGs (n = 5) (Figure 2 E). In 
individuals with three LAGs, the diameters of the outer, 
middle and inner LAGs, and absorbed zone ranged from 
282.7 um to 366.7 um, from 190.8 um to 269.4 um, from 
182.5 um to 214.4 um, and from 157.4 um to 197.7 um, 
respectively. The diameters of the outer and inner LAGs 
and absorved zone of the individual with two LAGs were 
288.8 um and 307.3 um, 173.0 um and 184.5 um, and 
145.0 um and 155.6 um, respectively. 


3.2 Application of skeletochronology to the Okinawaji- 
ma populations We constructed accurate methods to 
convert the number of LAGs into ages of individuals, 
based on collection date, SVL at capture, and horizontal 
diameters of LAGs (see Discussion). This method was 
applied for age determination to all geckos we handled. 

In the whole samples of G. kuroiwae, we found 0-3 
LAGs in juvenile, 3-7 LAGs in males, and 2-8 LAGs 
in females. Assuming that all individuals had hatched in 
September, the estimated maximal age was 70 months old 
in male and 83 months old in female. The longevity of 
G. kuroiwae was thus considered to be at least 83 months 
(ca. 7 years). 

With respect to the age structure of the population in 
the Motobu Peninsula, 28 out of 37 females (75.6 %) 
were younger than or as old as 48 months. The number 
of individuals dropped drastically in the age class of 
49-60 months, gradually decreasing thereafter (Figure 4 
A). The number of males decreased from 12 months, but 
this pattern may be merely due to the small sample size. 
The difference in age structure between the sexes was not 
statistically significant (median = 25.5 and 36.0 months 
in male and female, respectively; Mann-Whitney’s U test, 
Z = 1.327, P = 0.185). The scratch sample from scattered 
forests on the main body of Okinawajima mostly 
consisted of individuals younger than 48 months, except 
for one female of 53 months (Figure 4 B). 

Analysis using the von Bertalanffy growth curves 
demonstrated that asymptotic SVL was larger in female 
than in male in this species (male: asymptotic SVL + SE 
= 90.24 + 3.16 mm; female: 95.83 + 2.41 mm; Figure 5). 
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Figure 2 Cross-sections of second phalanx from the second digit of the left forearm. White and black arrows represent LAGs and hatching 
lines, respectively. White arrows with circles represent resorption lines. (A) A hatchling collected before the first winter (Group 1), (B—D) 
hatchlings collected after the first winter (Group 2), (E) a juvenile collected after the second winter (Group 3), and (F) a mature female esti- 
mated to be 83 months of age (Group 4). Right-bottom black bar indicates a length of 100 um for reference for all photos. 


On the other hand, no difference in growth coefficient 


was evident between the sexes (male: coefficient + SE 


= 0.047 4 


t 0.01; female: 0.047 4 


t 0.01). A statistically 


significant difference in the von Bertalanffy growth curve 
was observed between the sexes (F316 = 4.38, P < 0.01). 
Growth rates at 0, 12, 24, and 48 months were 1.61, 0.91, 


0.52, and 0.29 mm/month in male and 1.83, 1.05, 0.60, 
and 0.34 mm/month in female, respectively. 


4. Discussion 


4.1 Validity of skeletochronological age determination 
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We found a clear hematoxylinophilic line in hatchlings 
of Group 1 (Figure 2 A). Because this group included 
individuals collected from early to middle September, 
which appeared to be just after hatching, this 
hematoxylinophilic line is considered a hatching line 
(Guarino et al., 2010; Leclair and Leclair, 2011). 

Most individuals assigned to predefined age groups 2 
and 3 possessed two and three LAGs including a hatching 
line, respectively (Figure 2 B-C, E). Each of the LAG 
counts is considered as a complete number (= no loss) of 
LAGs in each group, because it is assumed that no LAGs 
are lost by endosteal resorption when the diameter of the 
absorbed zone is smaller than the minimum diameter of 
hatching line observed in Group 1 (165.0 um). Compared 
with Group 1, Group 2 possessed one additional LAG. 
This LAG would be formed during the first winter after 
an individual hatched (hereafter referred to as LAGI), and 
the inner LAG of this group corresponded to the hatching 
line (Figure 2 B-C). Then, compared with Group 2, 
Group 3 further possessed one additional LAG. Given that 
an additional LAG would be formed in the second winter 
(hereafter referred to as LAG2), no LAGs are expected 
to be formed during the summer months. The middle 
and inner LAGs in Group 3 corresponded to LAGI and 
hatching line, respectively (Figure 2 E). Consequently, 
it is considered that a single LAG is formed annually in 
winter. 

As we mentioned above, individuals with 
absorbed-zone diameter > 165.0 um may have 
lost LAGs through endosteal resorption. However, 
three out of five individuals in Group 2 and 3, 
even though the diameters of the absorbed zone 
were smaller than 165.0 um (145.0—-162.5 um), 
had LAGs one fewer than expected (Figure 2 C). 
The most plausible cause is a complete overlap between 
the hatching line and LAGI. Among individuals with 
both hatching line and LAGI in Group 2, thickness of 
the growth zone between the two lines varied (6.4-69.8 
um; Figure 2 B-D). These lines were sometimes partly 
overlapped, most likely due to insufficient growth 
of the individuals until the first winter (so called 
“rapprochement”; Castanet et al., 1992). An extreme case 
of rapprochement, in which LAGI was formed but just 
indistinguishable from hatching line, is most likely the 
case in the individuals that had fewer numbers of LAGs 
than expected. 

In the same manner, delayed or even no bone growth 
may have occurred in the single individual in Group 
2, in which no LAG was observed. Indeed, this was 
the smallest animal in the present sample (SVL = 37 


mm), similar in the size of a newly hatched individual 
(Tanaka and Nishihira, 1989), although it was caught in 
late February. Thus, although the hatching line/LAG1 
complex probably have been formed in this individual, 
it may be hard to recognize due to its location in the 
outermost layer of bone and direct contact with the well- 
stained membrane tissue. 

At the application of skeletochronology to routine 
age determination, an occasional overlap of LAGI with 
the hatching line as well as loss of LAGs by endosteal 
resorption will decrease the accuracy of this method. 
To address these problems, we attempted to distinguish 
hatching lines and LAGs based on the horizontal 
diameters of these lines. Figure 3 shows variations in 
horizontal diameters of LAG2, LAGI, hatching line, and 
absorbed zone. Diameter of absorbed zone was shown for 
all individuals (Groups 1—4), but diameters of hatching 
lines and LAGs were measured only for immature 
individuals with maximum LAG number in each age 
group 1-3. 
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Figure 3 Horizontal diameters of LAG2, LAG1, hatching line 
(HL), and absorbed zone by endosteal resorption (ER: resorption 
line or marrow cavity, whichever is larger). 


Diameters of the hatching line (range: 159.2—216.3 
um) and LAGI (range: 179.9-269.4 um) often 
overlapped, as shown in Figure 3. In contrast, LAG2 
(range: 282.7—-366.7 um) was obviously larger than LAG] 
and hatching line in all cases, and no overlapping was 
evident. Thus, LAG2 was distinguishable from LAG1 
and hatching line by its horizontal diameter. Horizontal 
diameter of the absorbed zone also overlapped with those 
of hatching line and LAGI, but never with that of LAG2. 
Based on these results, we formulated accurate methods 
of LAG counting and age determination as follows: (1) 
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Figure 4 Age structure of population in the Motobu Peninsula (A) 
and age composition in the pooled sample from discontinuous forests 
in the main body of Okinawajima Island (B). Gray, black, and white 
bars indicate juvenile, male, and female, respectively. 


when the diameters of all LAGs are < 270 um (including 
cases in which only one LAG is observed), LAGs include 
the hatching line, LAGI, and/or the hatching line/LAG1 
complex, and the age of individuals can be determined 
on the basis of the collection date and SVL at capture; (2) 
when the diameters of one or more LAGs exceed 280 um, 
the innermost one is LAG2, and the number of LAGs can 
be determined by counting outward from LAG?2. 


4.2 Longevity, age structure, and individual growth of 
the Okinawajima populations Goniurosaurus kuroiwae 
begins hatching from September and attains sexual 
maturity in the second winter after hatching (20 months), 
according to Tanaka and Nishihira (1989). In this study 
the potential longevity of this gecko was revealed to reach 
at least seven years. The longevity of gekkotan lizards 
estimated by various methods varies considerably. Short- 
lived species, such as Rhynchoedura ornata, live as long 
as 2 years (Read, 1999), while other species live much 
longer, such as Hoplodactylus stephensi (16 years; Hare 
and Cree, 2005), Homonota darwini (17 years; Piantoni 
et al., 2006), and Hoplodactylus maculatus (36 years; 
Bannock et al., 1999). The longevity of G. kuroiwae 
is modest among gekkotan lizards studied so far and 
comparable to those of other species distributed in the 
subtropical to tropical regions, such as Sphaerodactylus 
vincentii ronaldi (4 years; Leclair and Leclair, 2011), 
Teratoscincus roborowskii (> 6 years; Li et al., 2010), 
Oedura reticulata, and Gehyra variegata (> 9 years 
in both species; Kithener et al., 1988). Based on the 
analyses of age structure of two Okinawajima samples, 
most individuals appeared to live up to the age of 37—48 
months in Okinawajima (Figure 4), though this gecko 
potentially lives to approximately twice that age (83 
months). 
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Tanaka (1996) mentioned a low reproductive 
efficiency in G. kuroiwae, such as slow sexual maturity 
and small number of eggs produced. Moreover, our 
estimate of the age structure in this gecko suggests that 
it lives for only two or three breeding seasons. In the 
Motobu Peninsula, the situation looks better because 
females older than 48 months still accounted for 25% 
of females examined (Figure 4 A), and their life span 
therefore have encompassed four or more breeding 
seasons. This greater longevity, and consequently more 
chances of reproduction, is considered to contribute to 
the persistence of the population in this region. Expected 
positive correlation among age, female body size, and 
egg quality and/or clutch frequency may affect population 
persistence there, too (see below discussion). In contrast, 
the samples from the main body of Okinawajima included 
no older individuals of either sex (Figure 4 B), indicating 
a relatively low lifetime fecundity in this area. Due to this 
low fecundity as well as heavy habitat fragmentation, we 
suggest that the risk of local extinction is higher in this 
part of Okinawajima. 

With respect to individual growth of G. kuroiwae, the 
larger asymptotic SVL in female (90.2 mm in male vs. 
95.8 mm in female) was indicated by the von Bertalanffy 
growth curve (Figure 5), in agreement with a previous 
observation on female-biased sexual size dimorphism 
(Tanaka and Nishihira, 1987). The reproductive success 
of female G. kuroiwae, whose clutch size is fixed to 
two, may largely depend on the quality of each egg 
and frequency of clutch. In many lizards, egg size is 
positively correlated with the body size of offspring, and 
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Figure 5 The von Bertalanffy growth curves of male (solid line) 
and female (broken line) G. kuroiwae. Crosses, solid circles, and 
open circles represent juveniles, males, and females, respectively. 
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subsequently, with the performance of offspring (e.g., 
Sinervo, 1990; Warner and Andrews, 2002; Warner and 
Shine, 2009). Additionally, positive correlation between 
body size of female and egg size (mass) was reported in 
eublepharid geckos at both interspecific and intraspecific 
levels (Kratochvil and Frynta, 2002). Thus, larger 
body size is considered to benefit production of high 
performance offspring, increasing fitness of females of G. 
kuroiwae. 

In general, differences in SVL between male and 
female is caused by differences in age structure, growth 
rate, and/or timing of growth deceleration (Smirina and 
Tsellarius, 1996; Guarino, 2010; Leclair and Leclair, 
2011). Our results demonstrated no significant differences 
in age structure (25.5 and 36.0 months in male and 
female, respectively) and the von Bertalanffy growth 
coefficient (0.47 in both male and female), but showed 
difference in growth rate. Therefore, the faster growth 
rate in females is considered as a major cause of female- 
biased sexual size dimorphism in G. kuroiwae. Because 
most females were considered to die by 48 months of age 
and can experience breeding seasons only twice or thrice 
(Figure 4), females that attain larger body size until the 
first a few reproductive seasons are expected to achieve 
greater success. This may drive the faster growth rate in 
female. 
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